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Abstract 


Starting from the Braginskii equations, relevant for the tokamak edge region, a complete set 
of nonlinear equations for the geodesic acoustic modes (GAM) has been derived which includes 
collisionality, plasma beta and external sources of particle, momentum and heat. Local linear 
analysis shows that the GAM frequency increases with collisionality at low radial wave number kr 
and decreases at high kr- GAM frequency also decreases with plasma beta. Radial profiles of GAM 
frequency for two Tore Supra shots, which were part of a collisionality scan, are compared with 
these calculations. Discrepency between experiment and theory is observed, which seems to be 
explained by a finite kr for the GAM when flux surface averaged density (n) and temperature (T) 
are assumed to vanish. It is shown that this agreement is incidental and self-consistent inclusion 
of (n) and (T) responses enhances the disagreement more with kr at high kr ■ So the discrepancy 
between the linear GAM calculation, (which persist also for more “complete” linear models such 
as gyrokinetics) can probably not be resolved by simply adding a finite kr- 
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I. INTRODUCTION 


Common wisdom in fusion plasma science is that the transport of heat and particles 
in tokamaks are largely due to micro-turbulence driven by background gradients of density, 
temperature, momentum etc. The turbulence saturates via mode coupling, and in particular 


by interactions with self-generated large scale flow structures such as zonal flows[l|,and in 
some cases, especially near the edge, with geodesic acoustic modes (GAMs) 2-^. GAMs are 
an important class of oscillating zonal flows that appear due to toroidal geometry (i.e. due 
to geodesic curvature), and are easily observable in tokamak experiments due to their hnite 
frequency. They are usually classihed as an m = n = 0 perturbation in potential coupled 
with an m = 1, u = 0 perturbation in density or pressure, where m and n are the poloidal 



hnite frequency (usually a few kHz), distinct from that of broadband turbulence, GAMs are 
easier to detect, and thus have been observed on several tokamaks such as ASDEX Upgrade 
(AUG) [2l| using Doppler backscattering (DBS), TEXTOR [2^using 0-mode correlation 


rehectometer, and DIIID 


23| using beam emission spectroscopy (BES). As of today, GAMS 


are observed in the majority of the tokamaks in the world including recent observations of 
GAMs in Tore Supra 2^ using a DBS system. The common aspect of these measurements 
is that the GAMs are most prominent in the edge region, right inside the last closed hux 
surface, and extend into the near edge region (sometimes called the no man’s land due 


to a seemingly systematic discrepency between simulation and experiment 


25|). While 


gyrokinetics is accepted widely as the most general formulation for strongly magnetised 
plasmas of tokamak fusion devices, the applicability of gyrokinetic vs. huids models is still 
somewhat open to debate in the edge region. For most existing tokamaks, as one goes from 
the core to the edge, the collisions start to play a role, and the parallel connection length 
increases (since the safety factor q increases), dissipative drift waves, or resistive balloning 
modes, start to become important, therefore the validity of a huid description including the 
effects of collisions may actually be justi£ed|2a, l27 |. 

In this spirit, here we will dev elop a simple two fluid model for the description of the 


GAM, using Braginskii equations 


28 


29j within a drift expansion, in order to include the 


2 



























effects of collisions. In particular we include equations of contunity, momentum and heat 
for ions and electrons (i.e. using a generalized Ohm’s law for electrons) coupled with the 
Ampere’s law. The formulation allows us to include n||, Ay, Tj and perturbations of the 
GAM in a full set of nonlinear equations, which can be linearized and solved to obtain the 
GAM frequency including the effects of collisions and hnite /3, and hnite radial mode number 


The computed frequency is then compared with the radial prohle of GAM frequency 
that is observed in Tore Supra during a collisionality scan (assuming ~ 0). There is an 
apparent, systematic discrepancy between the theory and the experiment, which seems to 
be explained when a hnite kr is introduced for the GAM calculation assuming hux surface 
averaged density (n), ion temperature (Tj) and hux surface averaged electron temperature 
(Te) to be zero. However we believe that this agreement is incidental since it breaks down 


when higher harmonics (m = 2 etc.) are included in the calculation 30|, ImI]- The agreement 


also breaks down on self consistent inclusion of (n), (Ti^e) responses on the GAM dispersion 
for m = 1. This indicates that the discrepancy between the linear GAM calculation and 
the experiment, which persist also for more “complete” linear models such as gyrokinetics, 
is probably signihcant and can not be resolved by simply adding a hnite kr- 

The remainder of this paper is organized as follows. The complete set of nonlinear elec¬ 
tromagnetic equations with collisionality are obtained in SectioiilT] from the drift reduction 
of Braginskii equations. The fully nonlinear equations for GAMs are obtained in Section ITTn 
by taking appropriate hux surface averagings of the drift reduced electron and ion equa¬ 
tions. Linear GAM dispersion properties are obtained in Sectioi jTlI Al and comparison with 
experimental data are presented in Sectio JlIBI Finally the paper is concluded in Sectioi ilVl 


II. NONLINEAR MODEL EQUATIONS 


In order to formulate the nonlinear theory of electromagnetic geo desic acoustic modes 
(GAMs), we start with the simple two huid Braginskii equations 2^, [2^, where we keep 
the following: (i) the non adiabatic electron response with 6Te ^ 0 , and the electron- 
ion collisionality Vei- The model equations for GAMs are then derived from the density, 
momentum and temperature equations for each species j(= i, e). 

On 


dt 


+ V • (njVj) = 0 


( 1 ) 
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mini 


( 2 ) 


dvj 

dt 


+ Ui • V Ui = -Vpi - V ■ Hi + e \ E + -Vi X B \ + R 


0 = -Vpe - e X + Rei (3) 

Tj + Pj V • = -V ■ qj (4) 

where 

3 tiT''' —* 

Qj = 0.71nTjU\\ - K||V||Tj + K±'V±Tj + q^j + -z/j— -b x (5) 

2 cUcj 

Here the mass of the electron is neglected and E = V0 — {l/c)dA/dt. The collisional 
momentum transfer term is given by Rei = —Rie = neri^Jw — 0.71nV||Te, where J\\ = 
en{vi—v\\e), and U is the relative velocity between species j and i. The thermal conductivities 
for electrons are given by K||e = 3.16neTe/meVe, ^±e = 4.66neTeZ/e/mea;^e ^^cl for ions K\\i = 
^.QniTi/miUi, K,±i = 2niTih'i/miU)7. The diamagnetic heat flux is taken as q^j = b x 

VT,. 

In order to develop a drift expansion, we consider ion and electron perpendicular drift 
velocities in the low frequency regime {oj « ojd, « Ud] cv, oJd = eB/miC mode 
frequency, ion cyclotron frequency, respectively). These drift velocities consist of thei7 x B 
drift, the ion and the electron diamagnetic drifts, the ion polarization drift: 


ve = {c/B^)B X 
77*pi = {c/eniB^)B x V5pi 
v*pe = —{c/eneB^)B x V5pe 

Vp. = -£- (I + (-i. + ■ v) 


( 6 ) 

(7) 

( 8 ) 
(9) 


For the equilibrium scale lengths that are larger than the perturbation scales {i.e.,k±L < 1), 
we can separate the equilibrium (/o) and the fluctuating parts (Sf) in the above set 
of equations as / = /o + Sf. The complete set of resulting reduced nonlinear equa¬ 
tions for the perturbations {6f) are provided in the next subsection which is written 
in the following normalization scheme. The space time scales are normalized as r = 
r/ps, V|| = T„V||, t = tCs/Ln- The held quantities are normalized to their mixing 
length levels: (j) = {e5(j)/Te){Ln/ Ps), rii = {5ni/no){Ln/ ps), v = {5v\ii/cs){Ln/ Ps), Pi = 
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{6pi/Peo){Ln/Ps), A\\ = {2LnCs/(3psc){e6A\\/Teo). The remaining dimensionless parameters 
are : pi = L^JLrio, ^ + hO) T = Tio/Teo, /3 = SnPoe/B^, Lf = -dlnf/dx, 

Pe = -^rieoZ-^TeO) ^ = 0.51mer'ei-^n/''^iC5. The electron ion collision frequency is calculated 
from Uei = nZ'^lnX/{1.09xlO^^Te^^) where InX = 15.2—log{n/10^^)+log{Te) 32|. ps = CajoJci 
is the ion sound radius. The nonlinearities in the following equations originate mainly from 
the E X B drift nonlinearity i.e., vexb ■ V/ = [4>, f], the polariztion drift nonlinearity 
vexb ■ VV^^/ = [(j), 'V]_f] and the nonlinearity due to the parallel gradients with fluctuating 
magnetic helds, from V|| = V|| + 6Bj_ ■ V = Vfl — {(3/2) [y4||, ], Vj| being derivative along 

the equilibrium magnetic held. The (3 effects enter via perpendicular magnetic held line 
bending ehect through the expressions for and the instantaneous parallel derivative. 

Such models has also been used for edge turbulence simulations in the references 33l-l37(| 


A. Electron response 

When the drift expansion is considered in toroidal geometry, the electron continuity 
equation for density perturbation takes the form: 


1 f)(i) 1 f) 

M -rie- T,) - V,, (J|| - U||) 

= [0,Ue] - ^ [A||, J||-u||] (10) 


These and the following equations has been derived assuming large aspect ratio curcular 
hux surfaces. The second and third terms in the above equation results from the E x 
B convection of equilibrium density, the sum of divergence ol E x B diamagnetic drifts 
due to inhomogenous magnetic helds of the tokamak, respectively. The hrst term on the 
right hand is the E x B convective nonlinearity and the second term results from parallel 
derivative nonlinearity due to perpendicular magnetic huctuations. In the limit u << k\\Cs 
the perturbed parallel momentum equation for electrons reads: 


eJ\\ = -V|| {(j)-ne- 1.71Te) - ^ 


d . ^ 3 d 

- + {l + 1.71p,)-- 


^11 + f [" 4 ||, 0 -He- 1.71X1}) 
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Similarly, the electron temperature perturbation equation, in the same expansion, becomes: 




iTe 


) ^ 1 (9 d 

\n^) + -e„(cos 61-— + sin 6>—)Te + 
) 3 r od or 


2\ Idcf) 

3) r do 


1.07 


V T, 



+ 


1.07 


\ ^2r 89 


/3 dA\\ 




f bii.viiT. 


( 12 ) 


The hrst term on the right hand side is the E x B convective nonlinearity and the second 
term results from (5i?_Lpart of the parallel derivative. This equation has been derived using 
the continuity equation for V • v and the diamagnetic heat flux cancelation. 

Finally, the J\\ is related to the parallel vector potential via the Amperes law: 


Jii = -v5./i|| 


(13) 


B. Ion response 


The equations for the ion dynamics can similarly be obtained from the two fluid Braginskii 
equations using the drift expansion, with the aforementioned assumptions. The resulting 
equation for the continuity of ions take the form: 


drii 1 dcj) 
dt ^ r 80 


n 


A 8 . n ^ 

r 80 8 r 


(0 + TiTli + TiTi) 



[0, Hi] + V 


0 + Pi, V_L0 






vi0 


(14) 


Similar to the electron contunity equation, the second and third terms correspond to the 
E X B convection of the bacground density and the effects of inhomogenous magnetic held 
respectively, whilethe fourth term comes from the divergence of the polarization drift. The 
hrst term on the right hand side is the E x B convective nonlinearity, the second term is 
the polarization nonlinearity, and the third term results from parallel derivative nonlinearity 
due to perpendicular magnetic huctuations. Adding the electron momentum equation to the 
ion momentum equation and then using the electron temperature equation, one obtains the 
parallel ion velocity perturbation equation: 

18 8 \ 

cos6*-— + sin0—1 v\\ + V|| [r^ni + + Ue + T^] 

(t( 1 +0i) + 1+Ve) = - [4>,V\\] + ^ [A\\,Pi+Pe\ +niV|| {Pi+Pe) (15) 


8t 


^YjEr, 
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Here, the first term on the right hand side is the E x B convective nonlinearity, the second 
term results from perpendicular magnetic fluctuation induced parallel derivative nonlinear¬ 
ity, and the third term is the parallel acceleration term. 

Using the ion continuity equation for V • U and the ion diamagnetic heat flux cancelation 
the ion temperature perturbation equation becomes: 


a 2 
FtV‘-r' 


5 

-Er 

3 


Id . d \ ^ ( 2\ Idcj) 

+ (hi-gj -^ = 


i, Ti - -Hi 


(16) 


r d6 dr ^ 

Adding the electron and ion continuity equations after multiplying by respective charges 
and then assuming quasineutrality for the perturbations results in the plasma vorticity 
equation: 


d 


= -V 


4> + Ph V±0 


+1 

2 


^l|, ViA|| 


(18) 


where the f3 dependent nonlinear term comes from the perpendicular magnetic perturbations. 

The set of equations presented above for ions and electrons, provide a full drift-Braginskii 
system that can be used to describe the GAM oscillations including corrections due to finite 
/5 and collisionality, which may be relevant for the edge and near edge regions of tokamaks 
where GAMs have traditionally been observed. 


III. GEODESIC ACOUSTIC MODE 

GAMs are low poloidal mode number (m) axisymmetric fluctuations , that are supported 
by the geodesic component of the equilibrium magnetic curvature in tokamaks. In order to 
derive the set of equations that can be used to describe them, we start by taking the flux 
surface average of the vorticity equation flTT)) : 




(0) En-^ [(1 + Ti) {n sin 6) + n {Ti sin 6) + (Te sin 6)] 


dr 

0 + Pi, V±0 


+ 


/S 




(19) 


This shows that the flux surface averaged potential is linearly coupled to the m = 1 of density 
and temperature perturbations in the form of flux surface averaged {n sin 9) , (Tj sin 9) and 
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(Te sin 9) perturbations. This coupling happens due to the geodesic curvature. The effect 
of normal curvature for m = 1 mode is of the order of ps/r ~ ps/a ~ 10“^ and hence 
can be neglected. The nonlinear terms on the right hand side constitute the flux surface 
averaged poloidal momentum ffux/Reynolds stres and Maxwell stress and act as turbulent 
source/sink of vorticity. Equation (IT^ should be supplemented by the equation for (nsin^), 
which can be obtained by multiplying the ion continiuity equation by sin 6 followed by flux 
surface averaging: 


^ (nsin^) - (^sin^) - y Vr ((0) + n (n) + n (Tj)) - cos9) 

/S 


dt 


i, n] sin6') - ■ (j) + Pi,V±(j) sin + — ([A||, n||] sin 0) + (S'„ sin 6^) (20) 


Here, the term representing E x B convection of equilibrium density gradient is again 
dropped due to the fact that it is of the order of ps/r ~ p^/a ~ 10“^. Equation fl2U]) shows 
a linear coupling, this time with m = 1 of the parallel velocity ffuctution in the form of the 
flux surface averaged quantity (n|| cos6') . As is usually the case for GAMs, it is assumed that 
((0, n, Tj) sin^ 0) = ((0), (n), (Tj))(sin^ 6^), that is couplings to m = 2 and higher harmonics 
are ignored. In fact, the retention of m = 2 demands for the equations for m = 3 and so 
on, which continues up to infinity. This closure problem, and its possible resolutionwill be 
discussed in a future publication}^. The nonlinear terms on the right hand side acting 
as source/sink are poloidally asymmetric turbulent particle flux, poloidal momentum flux 
(i.e. Reynolds stress), and the electromagnetic component of parallel momentum flux. An 
asymmetry in the external particle source may also act as a source for the {n sin 9) component 
and therefore the GAM. The equation for (^v\\cos9'^ can be obtained by multiplying the 
parallel ion velocity equation (1T5|) by cos 9 followed by flux surface averaging: 

— (t>|| cos 9) + [(1 + a) + r* (Tjsin^) + (TeSin6')] = -([</>, t>||] cos 6*) 

([A||,pi +Pe] cos6>) + ((nV|| (pi +Pe)) cos9) + {SyCos9) (21) 

where the normal curvature term and second harmonic terms like (n|| sin 29) are again 
dropped. The various nonlinear terms on right hand side of the above equation can be 
identified as follows: The first term coming from the E x B convective nonlinearity is the 
COS 9 weighted, flux surface averaged, divergence of the parallel velocity flux. The second 
term is the electromagnetic analog due to perpendcular magnetic perturbation. The third 




term is the flux surface average of the turbulent parallel acceleration weighted by cos 6. 


This term survives only when there is a /cysymmetry breaking mechanism present which 
breaks the dipolar structure of acceleration in 6. The last term is the 6 symmetric part 
of the external velocity/ momentum source. The equation for (Tjsin^) , which appears in 
equations (IT^ and (1^ can be obtained by multiplying the ion temperature perturbation 
equation by sin 9 followed by flux surface averaging: 


dt 


{Ti sin 9) — -{n sin 9) 
o 


5 En 
3 ^ 


■V. (T,) = 


^,Ti 


-n 


sin 0 \ + {St sin 9) (22) 


The first nonlinear term on the right hand side is the divergence of flux surface average of 
the sin0 weighted heat flux minus 2/3 times particle flux. The second term is the poloidally 
asymmetric part of the external heating. 

The electron temperature equation for (Te sin 9) reads: 


d 


dt 


— {Tesm9) - - {nsm9) + -—Vr{Te) + 


5 Er 


1.07 (e 


/ Ti - n 

’ 3 


sin6' 


3 2 
{St, sin0) 


(T.sm»>-0,7l|^v;y/l||cos»> 

(23) 


The first nonlinear term on the right hand side is the divergence of heat flux minus particle 
flux weighted by sin 9. The second term is the poloidally asymmetric part of external heating 
Multiplying parallel electron velocity equation by cos9 and then taking the flux surface 
average gives: 

2q 


vVi- 


2 dt 


j (^ii cos6') — ^ [((/sin6') — {nsm.9) — 1.71 (TeSin6*)] 


/S 


= -- ([A||,(/ - n - 1.71Te] cos6*) 


(24) 


The electromagnetic nonlinear term on the right hand side comes from the perpendicular 
magnetic perturbation from the parallel derivative. Finally the equation for (0sin0) is 
obtained by multiplying the vorticity equation by sin6' and then taking the flux surface 
average: 

-- , 2 


(0sin6') - ( 1^11 cos6>) + [(1 + n) (n) + r* {Ti) + r* (Te)] 


= -(V- 


+ V_L(/ 




+ - {[AiiylAii] sm9) 


(25) 


The above equations are complemented by the transport like equations for (ue) and {Ty 
[(0sin0) - (uesin^) - (TeSin^)] = ([0,ne]) - ^([7l||, J|| -n||]) + {Snl26) 
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{TiSinO) = 


Ti - -n 


+ {^Ti} 


(27) 


^((^e) - ^ (n)) + ^v,. (TeSin0) = - 


- -n 


+ {Sn) 


(28) 


The 


and 


differ from the usual transport equations by the presence 


of a curvature term. This is due to the fact that these equations are obtained on taking 
flux surface average of the drift reduced equations fllOlllbp aud flT^ respectively where as the 
regular transport equations are obtained by taking the flux surface average of the starting 
fluid equations ([IM]). As a consequence (ue) and (Tj^g) also vary in GAM time scale rather 
than on slower transport time scale. The electron density and temperature equations arise 
because of the hnite beta extension which demands non-adiabatic electron response. In 
electrostatic case with adiabatic electron response one still needs to keep the equation for 
{Ti) to be consistent. Many of the previous papers except Ref. [^| on GAM happened to 
miss this somehow. Retention of this equation leads to another lower frequency branch with 
frequency going to zero when —)■ 0 which is distinctly different from the standard GAM 
whose frequency remains non-zero when —)■ 0. 

In the above equations S'„, St and Sy represents particle, heat and momentum sources 
respectively. This system of equations can be used to describe the complete nonlinear 
dynamics of GAMs using a reduced drift-Braginskii description including the effects of hnite 
/3 and collisionality. It is evident from the above equations that the GAM does contain 
m = 1 electromagnetic component which is in contrast to the previous Refs 


in line with the Refs 


40 


42 


4l| but 


43|. However these calculations can not explain the dominance 


of m = 2 electromagnetic perturbation as observed in some experiment and simulations 


44 


45| as the above equations are terminated at m = 1. Extension to m = 2 and beyond 


is deffered for future. The linearized form of the above equations can be written as a linear 
vector equation for the GAM state vector G 




(29) 


where M is a coupling matrix with elements depending on kr,Ti, and q. The matrix M is 
provided in the appendix by equation flAlj) . Upto m = 1 the GAM state vector G is made 
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of 


G = ((n), (0), (Tj), (Te), (rising), (^sin^), {A^ cos 9), (u|| cos 9), {Tism9), {TeSm9)) (30) 

and MisalOxlO matrix. In the z/ —)■ 0 bnt /3 7 ^ 0 limit it is straightforward to see that 
(TgSin^) = 0 and hence (Tg) = (n). Hence the GAM state vector will become 

G = {{n), (0), (Tj), {nsm9), {(psm9), (A|| cos9), (n|| cos9), {Tism9)) (31) 

and M becomes a 8 x 8 matrix. In the limit —)■ 0 and /9 ^ 0 from the above equations it 
follows that {nsin9) = (0sin0), (A||cos0) = (TeSin 6 ') = (n) = (Tg) = 0. Hence the GAM 
state vector reduces to 


G = (( 0 ), (Ti), (n sin 6 ^), (n| I cos 6 ^), (T^ sin 6 ^)) 
and M becomes a 5 x 5 matrix. 


(32) 


A. Linear GAM dispersion 

Linearizing the above set of equations fll9tl25p . taking Fourier transforms and neglecting 
(n), (Tj^e) , one obtains the following dispersion relation 


(jj 


1 - 


2 


_ kl{en/2qf (2.14cu + f(l.Q7/z/)(g^/2g)^)_ 

{{l3u/2 + iuk^oj - {en/2qY) {u + i(1.07/z/)(£„/2g)2) - l.2lklu{enl2qf 


Yq 


1 + 


5rj 


+ 


(2ci;/3) ((/9a;/2 + ii>kl)oj - (£„/2g)^) + Q.71klu{en/2q)‘ 


((/5a;/2 + wk‘^)u: - {en/2qY) [u + i(1.07/z/)(£„/2g)2) - l.2lklu:{en/2qY 
In the limit z/ — )■ 0 and /9 ^ 0 the above GAM dispersion relation becomes: 

-2 /, \ ^ \ / 


u‘ 


(i + G = R + 


2q 


1 + 


(33) 


(34) 


consistent with the basic GAM frequency as obtained by many authors (e.g. Ref 20|, l39|). 
However it is slightly different in the temperature ratio dependence as obtained fom the 
gyrokinetic calculations 


46 


48[ due to anisotropic temperature perturbations. 


In contrast fl33|) includes the effects of hnite jS and hnite collisionality. Notice that these 
effects appear together with kr in the above dispersion relation, implying that for k^ = 0 , 
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Figure 1. GAM frequencies vs kr with /3 and v as parameters. Color coding: blue(/3 = 0) to 
red(/? = 0.04) in steps of 0.004. Line styles; solid for z/ = 0 and dashed for u = 0.1. t = 1, q = 4. 


they can actually be neglected. We now explore some of the charecteristics of the GAM 
frequency, especially its scalings with /9, u and kr via the numerical solution of the dispersion 
relation fl33|) . Fig{T] shows the dispersion properties of GAM. Without collisionalities the 

h 0 at any kr- 


frequency decreases with kr monotonically and frequency also decreases wi 


This behavior is also consistent with the gyrokinetic calculation in Ref. 


41 


43]. At hnite 


collisionality, u = 0.1, the GAM frequency shifts up. The amount of upshift depends on kr 
at any given (3. At low values of (3 the upshift in frequency is more towards —?■ 0 and 
kr ^ 1 then when kr —)■ 0.5. Whereas at higher (3 values (> 0.02) the upshift in GAM 
frequency is noticable only beyond kr = 0.4. 


The radial group velocity at any f3 and kr is either zero or negative when i/ = 0. Whereas 
at hnite u the radial group velocity may either be negative, zero or positive depending on 
the values of /3, z/ and kr- 


Self consistent inclusion of hux surface averaged temperature leads to the following mod- 
ihcation of GAM dispersion in collisionless electrostatic limit 


u/ (l + kl) 



(1 + Tj + TiR 


TsinOnsinO 


(35) 
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Figure 2. GAM frequencies vs kr with effect of (Tj) for = 0, /3 = 0 


where the response RTsinensine of (T sin 6*) to (nsin6*) 

B 2/3 

■^TsinOn sin 9 - i 7^2 

1 — 2-^TTsin0 

and the response Rrusino of (T) to (nsin0) is given by 

Rtu sin 9 RtT sin 9 Rt sin9n sin 9 


where Rtt sine represents the response of (T) to (T sin0) 


Rtt sine — — 


1 5en 
a; 3 


On switching off the response Rtt sine, the equation [35] reduces to the simpler equation [SH 
Solutions of equation fl35|) are compared with the solutions of equation flMj) . It is seen that 
with self consistent treatment of the flux surface averaged ion temperature equation the 
frequency becomes non-monotonous in k^] decreasing at small but increasing at large 
kr- This can yield both inward and outward group velocites depending on the valaues of 
kr- Moreover another lower frequency branch appears (red curve in Fig.([2|)) which increases 
with kj. with the vanishing frequency at kj. = 0. 

Inclusion of (T) and (n) responses in the more general dispersion relation equation [33] 
for hnite collisionality and beta is straightfoward but tedious. The most general dispersion 
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Figure 3. GAM frequencies vs kr with effect of (Tj) and (n) for = 0 and 0.1 for different values 
of /3 


relation in principle can be obtained from 


\XI-M\ = 0 


(36) 


where the imaginary part of the eigenvalue A gives the real frequency and the real part of A 
gives the growth rate. It can best be studied numericaly for the eigenvalues A of the matrix 
M. FiglS] shows the dispersion relation with the self-consistent treatment of the average 
density and temperature equations. It is seen that in general the frequency decreases with jS 
with the rate of decrese depending on the values of kj. and v. However for a given {3 value the 
frequency increases with z/ for 0 < /c,, < 0.58 and decreases with u for > 0.58. The major 
differece between the exact dispersion relations Figl3]and the approximate one in FiglHare 
the following. The collisionality may enhance or reduce GAM frequency depending on the 
value of kr in the exact case whereas collisionality has always up shifting effect on GAM 
frequency in the approximate model. Also in the exact model the down shifting effect of f3 
on GAM frequency is much less pronounced than in the approximate model. 
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B. Comparison with experiment 


We compare the theoretical GAM frequency as given by the dispersion relations 
and the with the experimental GAM frequency observed in the Tore Supra tokamak 
for two different values of collisionality (i.e. shots ^^45494 and 7^45511) Equilibrium 
prohles of density and temperature are shown in Fig HI The temperature ratio Tj is greater 
than 1 towards the edge but less than 1 towards the core. The density remains almost the 
same towards the edge in these two discharges. Radial prohles of collisionality and beta, 
which are calculated using these equilibrium prohles, are shown in FigJSl 


solid lines high coll. #45511, dashed lines low coll. #45494 



r/a 


Figure 4. Equilibrium profiles at two different collisionalities. Solid lines are from high collisionality 
shot 7^45511and dashed lines from low collisionality shot ^45494. 


The radial prohles of experimental GAM frequencies are compared against the theoretical 
values as obtained from equations (IMi) . (I34j) . (l36i) and from Sugama’s formula |46| for kr = Q 
are shown in Fig. [6l Notice that our theoretical frequency compares well with Sugama’s for 
both the shots but experimental frequencies are lower than both. The frequency increases 
inward from the edge, which is consistent with increase of temperature, but the absolute 
values of experimental frequencies are about 50% below the theoretical values, consistently. 
We discuss below if this is due to hnite kr- Moreover, the experimental frequency is higher 
in low collisionality shots than in high collisionality shots. This might give the impression 
that GAM frequency goes down with collisionality, which is opposite to the theoretical 
prediction shown in the Fig{T]? This is in fact probably due to the change in temperature 
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Figure 5. (a) Radial profiles of collisionalities, (b) Radial profiles of plasma beta. 


profiles between the two shots, rather than the effect of collisionality. 
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Figure 6. Comparisons between experiment and theory. 


The prohles of theoretical GAM frequencies with and without collisions differs only 
slightly towards the edge where collisionality is higher than in the core. This shows that the 
collisionality and beta values are too low to make an appreciable impact on the frequency. 
However the absolute values of GAM frequencies are higher for low collisionality shots than 
for high collisionality shots, which is consistent with the upshift in temperature prohles as 
seen in FigJU 
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Figure 7. Radial scan of GAM frequency with krPs as parameter for high coll, shots ^^45511 


If one tries to explain the observed discrepancy between the theoretical and the experi¬ 
mental frequencies by invoking a hnite kr for the GAM one notes that the GAM frequency, 
as predicted by Eqn[33]goes down with increasing kr at any radius as shown in FigsI7]and 
[HI However it goes up with kr at the edge since collisionality is felt stronger at shorter scale 
lengths. The frequency values computed by assuming a hnite kr intersect with experimental 
observations for different wave numbers at different radii. In order to match with experi¬ 
mental observations, the radial wavenumber of the GAM has to increase with radius. The 
overlapping region of kr for high collisionality shots is [0.7,1.5] and for low collisionality 
shots is [1.3,1.7]. 

However we think that such a prohle of kr is rather unrealistic and is probably not the 
explanation of the observed discrepancy. The reason being that on self-consistent inclusion of 
(n) and (Tj e) responses through equation fl36|) breaks the monotonically decreasing behavior 
of GAM frequency on kr- Another important reason for this conclusion is that, in fact when 
other harmonics of the GAMs are considered (i.e. m = 2, m = 3 etc.), which are linearly 
coupled to m = 1, the kr dependence may be observed to change substantially. The gap 
between experiment and theory increases with kr at high kr on self-consistent inclusion of 
(n) and (Tj^e) responses at any radius as can be seen in the hgures ([9]) and fITOl) . 
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Figure 8. Radial scan of GAM frequency with krPs as parameter for low coll, shots #45494 
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Figure 9. Radial scan of GAM frequency with krPs as parameter for high coll, shots #45511, with 
self-consist treatment of (Ti^e) and (n) 

IV. CONCLUSIONS 

Using the reduced Braginskii equations under the drift approximation, a set of nonlinear 
electromagnetic equations retaining plasma beta and electron ion collisionality were ob¬ 
tained. Appropriate flux surface averaging were applied on the resulting set of equations in 
order to derive a fully nonlinear set of equations for the GAMs. This approach clearly shows 
that the GAM perturbations consist of (0), (0sin6*), (? 7 ,sin 6 '), (Tj^eSin6'), (u|| cos 6 *),(A|| cos^) 
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Figure 10. Radial scan of GAM frequency with krPs as parameter for low coll, shots ^^45494, with 
self-consist treatment of (Tj^e) and (n) 


up to the first poloidal harmonic. These equations can be used for studying both the non¬ 
linear drive and the linear oscillation of the GAM in a collisional, electromagnetic model 


of the plasma edge. However, one needs a better closure 


50 


52| in order to describe the 


collisionless GAM damping, which can in principle be included in the current model. Note 
however that initial studies in Tore Supra suggest that GAM damping near the edge region 
is mainly collisional as well. 


Linearizing these equations, a general linear dispersion relation is obtained, which con¬ 
tains the effects of electron ion collisionality and plasma beta. Following results are obtained 
without (u) and (Tj^e)- At zero collisionality and beta, the GAM frequency monotonically 
decreases with kr- At hnite beta the GAM frequency shifts down preserving the mono¬ 
tonically decreasing nature with kr- At hnite collisionality and low beta the frequency is 
shifted up, decreasing in kr at low kr and increasing at high kr- The upshift in frequency 
with collisionality is more prominent at low and high kr at low beta but, at high beta the 
collisional upshift is seen only at high kr values. However with self-consistent responses of 
(n) and (Tj^g) the results are modihed to the following. At zero collisionality and beta, the 
GAM frequency decreases with kr at low kr and incerases at high kr- At hnite beta the GAM 
frequency shifts down preserving the same non-monotonic decreasing nature. However the 
amount of down shift is more prominent in a small wave number range centered around 
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kr ~ 0.4. At finite collisionality the frequency is shifted up at low kr and shifted down at 
high kr, still preserving the non-monotonic behavior of frequency decreasing in kr at low kr 
and increasing at high kr- The upshift in frequency with collisionality is more at low kr than 
the down shift at high kr- We argue that these linear trends for the GAM frequency should 
be retained in the nonlinear regime, since the nonlinearity acts mainly as a source term for 
the GAM, and does not modify its response. 


The theoretical GAM frequencies obtained this way, were then compared to two TORE 
SUPRA shots, which were part of a collisionality scan and a consistent discrepancy were 
observed where the absolute values of experimental frequencies were about 50% below the 
theoretical values for kr = 0. In fact, without the (n) and (Ti,e) responses the experimental 
GAM frequencies may be tailored to match the theoretical values by assuming hnite kr 
values in the range kr G [0.7,1.5] for the high collisionality shots and kr G [1.3,1.7] for 
the low collisionality shot. However since these values are rather large, and since the kr 
dependence is mainly a feature of taking only the hrst poloidal harmonic of GAM, it is argued 
that “hnite kr effects” probably does not explain the observed discrepancy. Also the self 
consistent accounting of (n) and (Tj,e) responses makes the experiment-theory disagreement 
even worst on increasing kr- 


Note also that while the experimental GAM frequencies are lower at high collisionality 
than at low collisionality, this observed trend is due to the change in the temperature prohles 
and not due to the change in collisionality and plasma beta. This implies that in order to 
scan the collisionality dependence of the GAM frequency, one has to keep the temperature 
constant. Since these shots were part of a collisionality scan for the conhnement, it was the 
other dimensionless variables, such as p* etc. who were kept constant. 


These results leave the question of the discrepancy between the GAM frequency mea¬ 
sured in Tore Supra and the theoretical predictions. Notice that the discrepancy is equally 
important if one uses a more complex gyrokinetic formula, which is shown in hgure |6l We 
think that future work should include higher GAM harmonics in a similar fluid model of the 
edge, which may hopefully resolve this discrepancy. 


20 


ACKNOWLEDGMENTS 


The authors are thankful to X Garbet, Y Sarazin, and G Dif-Paradalier, T.S. Hahm 
and P.H. Diamond for discussions. This work has been supported by funding from the 
’’Departement D’Enseignement-Recherche de Physique Ecole Polytechnique” and it ihas 
been carried out within the framework of the EUROfusion Gonsortium and has received 
funding from the European Union’s Horizon 2020 research and innovation programme under 
grant agreement number 633053. This work was also partly supported by the french “Agence 
National de la Recherche” contract ANR JGJG 0403 01. The views and opinions expressed 
herein do not necessarily reflect those of the European Gommission. The authors are also 
extremley thankful to the annonymous referee for careful refereing and the suggestions for 
the improvement of the paper. 


Appendix A: The matrix M 

The matrix M in the equation |29] is 


M = {Mn] M^-Mt-, Mn 

sin 6 ; sinOi Ma, cosOi cosOi 


IVIT, sin 9 : iviT^ 

sin0 J 


where 


AG (0, 0, 0, 0, 0; 0? O? 

= (0,0,0,0,zen(l + rj)//cr,0,0,0, Anp//CrGen//^r) 

Mt^ = (0, 0, 0, 0, —i2enkr/^, z2e„/Cj,/3, 0, 0, ibenK/S, —i2e„/Cr/3) 

AGe = (0, 0, 0, 0, —i2enkr/3, i2enkr/3, 0, 0, 0, —i'lenkr/3) 

AGsinfl ( '^(^nkr/2^ ifinkj. j2^ Q^ i6^fcj./2, 0, 0, /2(5', 6 jj/2^, 0, 0) 

Af^sin^ ( ^^n(l T '^i)/2<k‘i’^ 0, tdfiTi f 2k‘p ^ %(i'fil2ky‘^ 0, 0, (ln/2Q^ 0, 0, 0) 

AG|| cose = (0, 0, 0, 0, en/fiq, -en//3g, -i/fc^2//3, 0, 0,1.71e„//3g) 

Af^ii cose = (0, 0, 0, 0, -e„(l + Ti)/2q, 0, 0, 0, -e„ri/2g, -en/2q) 

sin 0 ( '^^nkr / 3 ^ 'Id'fik'p j 3,0, Xdjikfl /6, 0, 0, d^ky, 1.38/2^, 6^/6^, 0, 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

-1.07(ej2g)tM0) 
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